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The electrical resistivity p of a metal is usually in- 
terpreted in terms of the average distance I an electron 
travels before it is scattered. As the temperature T is 
raised, p is increased and the apparent mean free path 
I is correspondingly reduced. In this semi-classical pic- 
ture / cannot be (much) shorter than the distance d be- 
tween two atoms. This has been confirmed forrnKiany 
systems and was considered a universal behaviour .Ein Re- 
cently, some apparent exceptions were found, e.g., alkali- 
doped FuUerenes&IO and high temperature superconduc- 
tors. These systems could, however, be in exotic states 
where only a small fraction of the conduction electrons 
contribute to the conductivity, and I would then have 
to be correspondingly larger to explain the observed p. 
Here we report on an essentially exact calculation for a 
model of alkali-doped Fullerenes, where the electrons are 
scattered by intramolecular vibrations. The resistivity at 
large T corresponds to I d, demonstrating that there 
is no fundamental principle of the type I > d. At large T 
the semi-classical picture breaks down and the electrons 
cannot be described as quasi-particles. 

We have also calculated the resistivity due to electron- 
electron scattering for a half-filled Hubbard model. In 
this case the resistivity saturates and / is not very much 
smaller than d. This difference is traced to the differ- 
ence between bosons and fermions. The resistivity is of- 
ten calculated using the Boltzmann equation. Although 
this equation is usually derived semi-classically, assum- 
ing I ^ d, va. our model for electron-vibration scattering 
it does not break down qualitatively at large T, where 
I <C d. For small T the calculated p due to electron- 
vibration scattering has a linear dependence on T and 
a strong depjendence on the pressure, in agreement with 
experiment E 

For A3C60 (A= K, Rb) at T - 500 K, the resistivity is 
p ~ 2 — 5 mficm. Assuming a spheijical Fermi surface and 
three conduction electrons per site,El we find that Z '--^ 1 — 2 
A is almost an order of magnitude smaller than the sepa- 
ration cZ = 10 A of the Ceo molecules. The experimental 
data have substantial uncertainties, but this is unlikely to 
influence the qualitative discussion of I. Different experi- 
mental methods (direct and optical) for different types of 
samples (thin films and doped single crystals) all suggest 
that / <C c? for large T. It has been suggested that the 
large p is due to the electrons being scattered inside the 
molecules. This would imply scattering between bands 
separated by at least 1 eV, and it plays a fairly small 
role for experimental T. In addition, we find that this 
interband scatteringB reduces the resistivity by providing 
an additional channel for conduction. 

The conduction in A3 Ceo takes place in a partly filled 
tiu band. The T-dependent part of the resistivity is as- 



sumed to be due to scattering against phonons (vibra- 
tions) with Hg symmetry. We therefore consider a model 
with a three-fold degenerate t\u level and a five-fold de- 
generate Hg Jahn- Teller phonon on each molecule, the 
hopping between the molecules and the coupling between 
the eleclmns and the phonons. The hopping takes inta 
accountB'E£l the orientational disorder of the molecules.EJ 
The one-particle band width is W = 0.6 eV and the 
phonon frequency is ujph- The coupling is determined by 
the dimensionless A = {5/3)N{0)g'^ /ojph, where A^(0) is 
the density of states per spin, g is the electron-phonon 
coupling constant and ujph is the phonon frequency, j— . 

We perform a finite temperature QMC calculation,E3 
treating the phonons quantum mechanically. We calcu- 
late the current-current correlation function for imagi- 
nary times and make a transformation to real frequencies, 
using a maximum entropy method.tJ The QMC method 
has no "sign-problem" , and the resistivity of the model 
can be calculated essentially exactly down to quite small 
T. 
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FIG. 1. The resistivity as a function of T. Results 
of the full QMC calculation, the Boltzmann equation 
(Bloch-Griineisen) and the bubble diagram are shown. The 
symbol x shows the p{T — 0) due to the orientational disor- 
der. The figure illustrates that p can become extremely large, 
that the bubble calculation is quite accurate and that there is 
no qualitative break-down of the Boltzmann equation at high 
T. 

Fig. |l| shows the resistivity for a cluster of 48 Ceo 
molecules for A — 0.5 and uiph = 0.2 eV. The QMC cal- 
culation (full line) shows that the resistivity can become 
very large, corresponding to / ^ 0.7 A at T — 0.5 eV. 
By also considering an unrealisticly large T, we empha- 
sise the lack of a limitation of the type I > d. Then 
there cannot be any general principle requiring / > d, 
because such a principle would then also apply to the 
present model. Qualitatively similar result were obtained 
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by MilliS|-ei aZ.,0 using the dynamical mean- field theory 
(DMFT)li3 and assuming classical phonons. Their calcu- 
lation, however, does not prove that I <^ d is possible, 
since it involves approximations. The moderate differ- 
ences to our results are probably due to different models 
and their use of classical phonons. 

While the QMC calculation is essentially exact, it is 
hard to interpret the results. For this purpose we use a 
diagrammatic approach. The idea is to simplify the di- 
agrams as far as possible without loosing the qualitative 
agreement with the QMC results. These diagrams are 
then studied to obtain a qualitative understanding of the 
results. This approach (Kubo formalism) requires the 
calculation of a bubble diagram including vertex correc- 
tions (see Fig. ||a). We neglect the vertex and calculate 
the bubble using the electron Green's function from the 
QMC calculation. The resulting resistivity (dashed line 
in Fig. ^ is practically identical to the QMC result, jus- 
tifying the neglect of vertex corceptions for the present 
model. It was shown by Holsteinli3 that in the limit of a 
broad electronic band, all vertex corrections except lad- 
der diagrams can be neglected and that a Boltzmann 
equation can be derived. Holstein's derivation is not valid 
for the narrow band considered here, but our calculations 
show that his arguments are still rather accurate. For our 
model with a q-independent electron-phonon coupling, 
even the ladder diagrams can be neglected. Essentially 
following Holstein we obtain approximately a Boltzmann 
like conductivity 
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where N{lu) is the density of states, / is the Fermi func- 
tion, S(w) is the electron self-energy, jk is the current 
matrix element for a state with the label k and the en- 
ergy Ek- We interpret Im S as the inverse of the relax- 
ation time. For a large T, Im S becomes comparable to 
or larger than the one-particle band width. Since Im E, 
is related to the inverse life-time of the quasi-particles, 
this means that the life-times become so short that the 
concept of a quasi-particle breaks down at large T. 
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FIG. 2. Relevant diagrams. The current-current re- 
sponse function (a) and two approximations to the elec- 
tron self-energy ((b) and (c)) are shown. The full and 
dashed lines represent electron and phonon Green's functions. 
Self-consistent Green's functions are used in (a) but not in (b) 
or (c). 

The resistivity thus depends crucially on S. To under- 
stand its behavior, we consider the diagram in Fig. ^ 
calculated with bare Green's functions and for simplicity 
neglecting the orbital degeneracy 
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is the Rose occupation number. For large T, ub becomes 
large, leading to a large Im E, a small a and a large p. 
The Bose nature of the phonons is therefore of crucial 
importance for our result. 
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FIG. 3. The resistivity of the nondegenerate Hubbard 
model. Results are shown for different values of the ratio be- 
tween the Coulomb repulsion U and the band width W . The 
figure illustrates how the resistivity saturates for large T in 
the electron-electron scattering model. 

This is further illustrated by comparing with the re- 
sistivity due to the electron-electron scattering. Using 
the DMFT we have calculated the resistivity for a non- 
degenerate Hubbard model on a Bethe lattice with d = 1 
A. We focus on the half-filled case, which is relevant for 
A3 Geo, aud we do not consider the case of a doped Mott 
insulator .113 The impurity problem is solved with a QMG 
method. Fig. ^ shows piT) for different values of the on- 
site Coulomb interaction U . For U /W < 1 the system is 
a metal and p{T) grows with T, while for U /W > 1 it is 
an insulator and p{T) decreases with T. The important 
observation is that in the metallic case p{T) saturates at 
p ~ 0.4 mflcm, which corresponds to l/d ~ 1/3. Thus, 
in contrast to the electron-phonon scattering case, I is 
not very much smaller than d in the metallic state of this 
model and within DMFT. 

To understand this, we study the electron sclf-encrgy S 
to second order in U, since E determines p in the DMFT. 
For low T, there is little scattering due to the small 
phase space available, as controlled by the Fermi func- 
tions. As T increases, the available phase space grows 
and p increases. However, for large T, p essentially sat- 
urates, since the Fermi functions approach a constant 
value. This is in strong contrast to the Bose occupation 
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numbers (Eq. (||)), which keep increasing with T. The 
qualitative difFerence between the two models for large 
T can then be traced to the difFerence between fermions 
and bosons. 

We now address the validity of the Boltzmann equation 
for the case of the electron-phonon scattering in view of 
I <C d. We have calculated the resistivity using the Ziman 
version otthe Bloch-Griineisen solution of the Boltzmann 
equationlla and added the resistivity due to thCiJirienta- 
tional disorder as a T-independent contribution.Ej Using 
the plasma frequency Up = 1.2 eV, we obtain the dotted 
line in Fig. |l|. The Boltzmann result is larger than the 
QMC result for large T, but there is no qualitative break 
down of the Boltzmann equation, although I d and 
the quasi-particle concept is not applicable. The justifi- 
cation for the Boltzmann equation in this limit is not the 
semi-classical derivation, but the (approximate) deriva- 
tion from the full quantum mechanical Kubo formulation 
(Eq. (|l|)). The proper language is not in terms of a very 
short mean free path, but in terms of a very broad spec- 
tral function (Im S large), as discussed above. 

The QMC calculation gives an approximately linear T 
dependence. This agrees with the exped.mental result 
that p is linear down to about 100-200 K.El although the 
linearity is not as pronounced as for some high temper- 
ature superconductors. The result may seem surprising, 
since at small T the probability of exciting finite energy 
phonons is exponentially small as is the contribution to 
p. Calculating the bubble diagram in Fig. ||a using elec- 
tron Green's functions with the self-energy in Fig. 
indeed gives an exponential behaviour. The use of the 
QMC Green's function, however, gives a linear behavior, 
in agreement with the full QMC calculation. The reason 
is that the QMC Green's function also involves processes 
like in Fig. where a virtual phonon is created followed, 
by the decay of this phonon into an electron-hole pair.E2l 
The excitation energy of such a pair can be atbitrarily 
small, leading to a quadratic T-dependence for pB3 In our 
model this goes over in an approximately linear behavior 
already for very small T. 

Calculations for different values of A show a strong A- 
dependence for small T. This is in agreement with the. 
strong pressure dependence observed experimentally,^ 
since the application of pressure increases the band width 
and reduces N{0) and thereby A. For A ^ 1 there is a 
transition to an insulating state. 

Our use of a fairly realistic model with a three-fold 
degenerate electronic level and a five-fold degenerate 
phonon is not crucial for the qualitative behaviour of 
p{T), like I ^ d for large T. An essential feature is, 
however, the use of intramolecular phonons acting on a 
narrow band. The phonon induced shifts (~ T) of the tiu 
levels can be substantial compared with the small band 
width. It would be interesting to perform calculations 
for a model where the phonons couple to the hopping 
matrix elements and where interband transitions become 
important for large T. Such a model would be more rel- 
evant for, e.g., transition metal compounds, where the 



resistivity typically saturates when I ^ d. 
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